Necessary and sufficient condition for non-zero quantum discord 
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Quantum discord characterizes "non-classicality" of correlations in quantum mechanics. It has been proposed 
as the key resource present in certain quantum communication tasks and quantum computational models without 
containing much entanglement. We obtain a necessary and sufficient condition for the existence of non-zero 
quantum discord for any dimensional bipartite states. This condition is easily experimentally implementable. 
Based on this, we propose a geometrical way of quantifying quantum discord. For two qubits this results in a 
closed form of expression for discord. We apply our results to the model of deterministic quantum computation 
with one qubit (DQC1), showing that quantum discord is unlikely to be the reason behind its speedup. 
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Introduction. — Quantum states of a composite system can 
be divided into entangled and separable once. Entangled 
states display "nonlocal features" violating Bell's inequali- 
ties |H and are considered a necessary resource for quan- 
tum communication and pure quantum computation allowing 
computational speedup over the best classical algorithm f^j]. 
On the contrary, separable states are generally considered as 
purely classical, since they do not violate Bell's inequalities 
and can be prepared by local operations and classical commu- 
nication. However, it is valid to ask if highly mixed states, 
and in particular separable states, are completely useless from 
quantum information perspective. Recent investigations give 
compelling evidences that this is not the case. A highly mixed 
state in the DQC1 model y[ is believed to perform a task 
exponentially faster than any classical algorithm ("without 
containing much entanglement"). Furthermore, it has been 
shown that even some separable states contain nonclassical 
correlations JUS] and can create an advantage for computing 
and information processing tasks over their classical counter- 
parts 0-Gi]. 



The "non-classicality" of bipartite correlations is measured 
via quantum discord ||4|]- the discrepancy between quantum 
versions of two classically equivalent expressions for mu- 
tual information. Recently, it has been shown that almost all 
quantum states have non-vanishing discord 11211 . Quantum 
discord was proposed as a figure of merit for characterizing 
the nonclassical resources present in the DQC1 llOtl . It has 
been shown that initial zero-discord system-environment state 
is necessary and sufficient condition for completely-positive 
map evolution of the system when the environment is traced 
out Il3l. [l4tl . Furthermore, in Ref. 1 15] is demonstrated that 
if the state can be locally broadcasted than it has vanishing 
discord. 

Despite increasing evidences for relevance of quantum dis- 
cord in describing non-classical resources in information pro- 
cessing, there is no straightforward criterion to verify the 



presence of discord in a given quantum state. Its evalua- 
tion involves optimization procedure and analytical results 
are known only in a few cases iflrjll . In this Letter we de- 
rive the necessary and sufficient condition for non-vanishing 
quantum discord. The criterion is simple and also experimen- 
tally friendly, since it can be evaluated directly from a (sub)set 
of measurements standardly used for quantum state tomogra- 
phy. Based on this, we introduce the geometrical measure of 
discord and derive an explicit expression for the case of two 
qubits. Finally, we give arguments putting in question appro- 
priateness of quantum discord to describe the non-classical 
resource in DQC1 computational model. 

Quantum discord. — Correlations between two random 
variables of classical systems A and B are in information the- 
ory quantified by the mutual information I(A : B) = H(A) + 
H(B) - H(A,B). If A and B are classical systems, than H{.) 
stands for the Shannon entropy //(p) = - 2/ Pi 1°S Vu where 
p = (p\,p2, • ■ • ) is the probability distribution vector, while 
H(., .) is the Shannon entropy of the joint probability distribu- 
tion pij. For quantum systems A and B, function H(.) denotes 
the von Neumann entropy H(p) = -Trp log p where p is the 
density matrix. In the classical case, we can use the Bayes 
rule and find an equivalent expression for the mutual informa- 
tion I(A : B) = H(A) - H(A\B) where H(A\B) is the Shannon 
entropy of A conditioned on the measurement outcome on B. 
For quantum systems, this quantity is different from the first 
expression for the mutual information and the difference de- 
fines the quantum discord. 

Consider a quantum composite system defined by the 
Hilbert space 'Kab = "Ha ® 'Hb- Let dimensions of the lo- 
cal Hilbert spaces be dim'H A = dA and dim7Y B = ds, while 
d = dim'7-/ AB = c/a^b ■ Given a state p (density matrix) of a 
composite system, the total amount of correlations is quanti- 
fied by quantum mutual information [l7\\ : 

I(p) = H(p A ) + H(p B ) - H(p), (1) 

where H{p) is the von Neumann entropy and pa,b = TrB,A(p) 
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are reduced density matrices. A generalization of the classical 
conditional entropy is H(pb\a), where p B \A is the state of B 
given a measurement on A. By optimizing over all possible 
measurements in A, we define an alternative version of the 
mutual information 

Q A (p) = H(p B ) - min {i y ^ PkH(p m ), (2) 



where p%\k = T*A(Ek 8 lBP)/TrtE* 8 11 bp) is the state of B 
conditioned on outcome k in A and {Ek} represents the set of 
positive operator valued measure elements. The discrepancy 
between the two measures of information defines the quantum 
discord J3H: 



D A (p) = I{p) - Q A {p). 



(3) 



The discord is always nonnegative |4J] and reaches zero for 
the classically correlated states 0. Note that discord is not a 
symmetric quantity D A (p) + Dgip) and D A refers to the "left" 
discord, while Dg refers to the "right" discord. The state p for 
which D A {p) = Dr(p) — is completely classically correlated 
in a sense of 11181 Il9ll From now on, when we refer to the 
discord we mean the "left" discord D A . 

To give an example of a state with non-vanishing dis- 
cord consider the two-qubit separable state in which four 
nonorthogonal states of one qubit are correlated with four 
nonorthogonal states of the second qubit: 
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4 -(|0)(0| 8 l+X+l + |1><1| 8 l-X-l + l+X+l 8 UXH 
+I-X-I® 10X01). 



(4) 



Unlike the state above, one can show that the state p is of 
zero-discord if and only if there exist a von Neumann mea- 
surement {IXt = |^)(^|} such that I20I1 



^(rr. 8 i B )p(ru 8 1 B ) = p, 



(5) 



In other words the zero-discord state is of the form p = 
T*k Pk\4'k){4'k\ ®Pk where W^k)} is some orthonormal basis set, 
Pk are the quantum states in B and pk are non-negative num- 
bers such that Yik Pk — 1- 

An easily implementable necessary and sufficient condi- 
tion. — Let us choose basis sets in local Hilbert-Schmidt 
spaces of Hermitian operators, {A„} and {B m } where n = 
1 . . .d 2 K and m = 1 . . . c/ 2 . We decompose the state p of 
composite system into p = Y^nm r nmA n 8 B,„. The coeffi- 
cients r nm define d 2 K x d\ real matrix R which we call the 
correlation matrix. We can find its singular value decompo- 
sition (SVD), URW T = diagfci, c 2 , . . . ] where U and W are 
d\ x d 2 A and d\ x d^ orthogonal matrices, respectively, while 
diag[ci, C2, . . . ] is d\ X diagonal matrix. SVD defines new 
basis in local Hilbert-Schmidt spaces S„ = U nn 'A n i and 
Fm = Yum 1 W mm >B m >. The state p in the new basis is of the 
form p = 2n=i c nS n ® F„ where L = rank/? is the rank of 
correlation matrix R (the number of non-zero eigenvalues c„). 



The necessary and sufficient condition © becomes 
En=i c„(Yjk U kS „^k)®F n = 2«=i c n S„®F n and it is equivalent 
to the set of conditions: 



2^njfeS n II t = S„, n = l...L, 



(6) 



or equivalently [^n,^] = for all k, n. This means that the 
set of operators {S „) have common eigenbasis defined by the 
set of projectors {rU}. Therefore, the set {11^} exists if and 
only if: 



[S n ,S m ]=0, n,m=l...L. 



(7) 



In order to show zero discord we have to check at most 
L(L - l)/2 commutators, where L = rank/? < minjof^, d^}. 
Now, recall that the state of zero discord is of the form 
P = Tifti PkUt ® Pk< therefore is a sum of at most dA product 
operators. This bounds the rank of the correlation tensor to 
L < dA- Thus, the rank of the correlation tensor is the simple 
discord witness: If L > dA, the state has a non-zero discord. 

Correlation matrix can be obtained directly by simple mea- 
surements usually involved in quantum state tomography. 
However, the detection of non-zero discord does not neces- 
sarily require measurement of all (d A dB) 2 elements of the cor- 
relation matrix (full state tomography). It is sufficient that the 
experimentalist measures that many elements of the correla- 
tion matrix until he finds dx + 1 linearly independent rows (or 
columns) of the correlation matrix. 

Geometric measure of discord. — Evaluation of quantum 
discord given by equation (01 in general requires consider- 
able numerical minimization. Different measures of quan- 
tum discord l2lll and their extensions to multipartite sys- 
tems 1l9h have been proposed. However, analytical expres- 



sion are known only for certain classes of states lfl6ll . Here 
we propose a following geometric measure 



Df(p) = min^njlp -^|| 2 , 



(8) 



where Oq denotes the set of zero-discord states and \ \X- Y\\ 2 = 
Tr(X - Y) 2 is the square norm in the Hilbert-Schmidt space. 
We will show how to evaluate this quantity for an arbitrary 
two-qubit state. 

Two-qubit case. — Consider the case "Ha = *Hb = C 2 . We 
write a state p in Bloch representation: 



.3 3 3 

p = -(181 + V x;a-;8l + Vyill8crj+ V Tijcr^o-j), (9) 



where x-, = Trp(<T,- 8 11), y,- = Trp(l 8 cr,-) are components of 
the local Bloch vectors, 7y = Trp(<r, 8 cry) are components of 
the correlation tensor, and cr,, i e {1, 2, 3}, are the three Pauli 
matrices. To each state p we associate the triple {x,y, T). Now, 
we characterize the set Qo. A zero-discord state is of the form 
X = Pil^iX^il ®Pi + P2l^2><^2l ®P2, where fl^i), \f 2 )} is a 
single-qubit orthonormal basis, p\2 are 2 x 2 density matrices, 
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and pi2 are non-negative numbers such that p\ + pi = 1. We 
define t — p\ —p% and three vectors 



if Ml), 
Tr(pipi ± pipi)&. 



(10) 
(11) 



It can easily be shown that te and s + represent the local Bloch 
vectors of the first and second qubit, respectively, while the 
vector sL is directly related to the correlation tensor which 
is of the product form T = es[. Therefore, a state of zero- 
discord x nas Bloch representation^ = {te, s + ,efF_), where 
1 1^1 1 - 1> ll*±ll ^ 1 an( l t e [ — 1> !]■ The distance between states 
p and x is given by 

Hp-;tII 2 = IIpII 2 -2Tw + iwi 2 

= \a + m 2 + ii^i 2 + \\t\\ 2 ) 

- —(1 + txe + ys + + eTsJ) 



(12) 



\s-\\ 2 ), 



where ||r|| 2 = TrT J T. First, we optimize the distance over 
parameters s ± and t. The function of equation fTHi is convex 
and quadratic in its variables t, s ± . It is straightforward to see 
that its Hessian is a positive and non-singular matrix. There- 
fore the function has a unique global minimum. The minimum 
occurs when the derivative is zero: 

llp-^ll 2 



1 



dt 
Wp-xW 2 

ds+ 

Wp-xW 2 



= -(-xe + t) = 0, 

= = 0, 



= -(-r 1 

2 



'+/_) =0, 



(13) 
(14) 

(15) 



which gives the solution t - xe, s+ — y and s*- = 
T T e. Since the solution lies within the range of param- 
eters, \xe\, \\y\\, ||r T e|| < 1 it represents the global mini- 
mum. After substituting the solution we obtain \\p - x\\ 2 - 
\ (\\x\\ 2 + ll^ll 2 - eixf + TT T )e) which attains the minimum 
when e is an eigenvector of matrix K = icx 1 + TT T for the 
largest eigenvalue. Therefore, we have: 



Df(p) 



1 

4< 



\x)\ 2 + \\T\\ 2 -k max ), 



(16) 



where k max is the largest eigenvalue of matrix K = xx^ + TT T . 
Next, we apply our criterion to a class of states. 

States with maximally mixed marginals. — We consider an 
example of two qubit states with maximally mixed marginals. 
Such a state is locally equivalent (under some local unitary 
transformation U\ ® U2) to a state p(t) — (1 ® 1 + 2^ =1 f;0"j ® 
<Ti)/4, where t = {t\,t2,h). The state p(t) is physical if /^be- 
longs to the tetrahedron (Figure [TJ defined by the set of ver- 
tices (-1,-1,-1), (-1, 1, 1), (1,-1, 1) and (1, 1,-1), while is 
separable if t belongs to the octahedron defined by the set of 
vertices (+1,0,0), (0,+l,0) and (0,0, +1) JH. Simple cal- 
culation shows that D (2) (t) = \{t\ + t\ + t\ - max{f 2 , f 2 , t\}). 
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FIG. 1 . The set of two-qubit states with maximally mixed marginals 
(i.e. the reduced states of individual qubits are completely mixed). 
Physical states belong to the tetrahedron, among which separable 
ones are confined to the octahedron. The zero-discord states are la- 
beled by the red lines (it is therefore clear that almost all states have 
non-zero discord [12]). The states with maximal value of discord 
correspond to the vertices of the tetrahedron (the four Bell states). 
Among the set of separable states, those which maximize discord are 
the centers of octahedron facets (±1, ±1, ±l)/3 (black dots). 



The zero-discord states have at most one non-zero component 
of vector t (FigureQ] red lines). The function D^ 2 \i) reaches 

(2) 

its maximal value of Dy = 1 /2 at the vertices of tetrahedron 
which represent the four Bell states. Within the set of sep- 
arable states (octahedron) its maximal value of Z?^ 2) = 1/6 
is attained at the centers of octahedron facets (±1, ±1, ±l)/3. 
They represent the states 



1 1 

Phhh = 7(1 ® 1 + o ^(-1)"^ ® cr k ), 
*=1 



(17) 



where = +1, and can intuitively be understood as equal 
mixture of "maximally non-orthogonal" states. The states are 
symmetric under exchange of subsystems, thus they have the 
same value of "left" and "right" discord Da = Dg- 

DQC1 model. — In J3t], Knill and Laflamme introduced the 
model of mixed-state quantum computing which preforms the 
task of evaluating the normalized trace of a unitary matrix ef- 
ficiently. The corresponding quantum circuit is shown in Fig- 
ure [2] The input state is a highly mixed separable state and 
consists of a control qubit in the state ^(11 +acr 3 ), where a de- 
scribes the purity, and a collection of « qubits in the maximally 
mixed state ^ 1„, where 1„ is the n-qubit identity. The DQC1 
circuit consists of the Hadamard gate applied to the control 
qubit and a control n-qubit unitary gate U n . The output state 
is: 

P=^L(li®ln + a\lX0\®U n + a\0)(l\®Ul). (18) 
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FIG. 2. The quantum circuit for estimating the normalized trace of 
the unitary matrix U„ using the model of deterministic computing 
with one quantum bit (DQC1). H stands for the Hadamard gate. 
The control (top) qubit is measured in the cr l and cr 2 basis, and the 
expectation values give the real and imaginary part of normalized 
trace r = Trt/„/2" with the overhead scaling as \ floTl . 



We consider only the cases a + 0, otherwise the state at the 
output is completely mixed and therefore cannot accomplish 
the task. After measuring the control qubit at the output in the 
eigenbasis of <j\ and o~%, we retrieve the normalized trace of 
the unitary matrix r = Tri/„/2" with the polynomial overhead 
scaling I /a 2 fioll . 

The control qubit is completely separable from the rest of 
the qubits. The output state has vanishingly small entangle- 
ment across any bipartite split that groups the control qubit 
with some of the mixed qubits 0], However, there is strong 
evidence that DQC1 task cannot be preformed efficiently us- 
ing classical computation The question is what brings a 
"speed-up" in the considered task? The quantum discord was 
proposed as a figure of merit for characterizing the resources 
present in DQC1 model Holl , It has been shown that for al- 
most every unitary matrix U„ (random unitary) the discord in 
the output state (fl8l l is non-vanishing. Here we derive an ex- 
plicit condition for characterizing the correlations in the out- 
put state and show that the discord is unlikely to be the source 
of speedup. We re-write it into a form: 



1 



li ® t„ + acr\ 



U„ + ul 



+ aa"2 i 



2i 



(19) 

Now, we apply the condition Q. The operators cr\ and <x 2 do 
not commute, therefore, the state p is of the zero-discord if and 
only if the operators u "+ u - and V "^" are linearly dependent, 
or equivalently Jj\ — kU„. This is possible if and only if U n = 
e'^A, where A 2 = 1 is a binary observable. For such a unitary 
all the correlations at the output of DQC1 circuit are classical. 
However, it is very unlikely that the normalized trace of e'^A 
can be evaluated efficiently on a classical computer, since all 
it's eigenvectors can be arbitrarily complex (random states). 

We emphasize that our measure of discord is not mono- 
tonic under local operations. This, however, is not a shortcom- 
ing, as discord, unlike entanglement and mutual information, 
can in fact increase as well as decrease under local operations 
(even without the presence of classical correlations). A sim- 
ple example of the local increase is to start from a zero-discord 



state |00)(00| + |11)(11| and transform, say the first qubit, so 
that |0) — » |^o ) and |1) — » \if/i), such that and are 
not orthogonal. The resulting state, lO^oXO^ol + 
clearly has a non-vanishing discord. Finally, we point out that 
our method can be extended to any number of subsystems, 
though evaluating the measure of discord becomes progres- 
sively more difficult with increasing number of subsystems 
and their dimensionality. 
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Note added in proof. — A related work was done by 
Datta H. 
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